Fluctuations in protein synthesis from a single RNA template: 
stochastic kinetics of ribosomes 
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Proteins are polymerized by cyclic machines called ribosome which use their messenger RNA 
(mRNA) track also as the corresponding template and the process is called translation. We explore, 
in depth and detail, the stochastic nature of the translation. We compute various distributions 
associated with the translation process; one of them, namely dwell time distribution, has been 
measured in recent single ribosome experiments (Wen et al. Nature 452, 598 (2008)). The form 
of this distribution predicted by our theory is consistent with that extracted from the experimental 
data. For our quantitative calculations, we use a model that captures both the mechano-chemistry 
of each individual ribosome as well as their steric interactions. We also demonstrate the effects of 
the sequence inhomogeneities of real genes on the fluctuations and noise in translation. In principle, 
our new predictions can be tested by carrying out in-vitro experiments. 

PACS numbers: 87.16.af; 87.16. dj 



A genetic message, chemically encoded in the DNA, 
is first transcribed into a messenger RNA (mRNA) from 
which it is then translated into proteins [l[ . Both mRNA 
and proteins are linear polymers of monomeric subunits 
called nucleotide and amino acid, respectively. The ge- 
netic code contained in the sequence of codons (triplets 
of nucleotides) on an mRNA is translated into the corre- 
sponding sequence of amino acids by a macromolecular 
machine, called ribosome [|J ■ A ribosome is a cyclic ma- 
chine. Each mechano-chemical cycle of this machine con- 
sists of several steps which result in the translocation of 
the ribosome by one codon on the mRNA template and 
the elongation of the protein by one amino acid. Thus, 
the mRNA template also serves as the track for motor- 
like movement of the ribosome during translation 

Si- 

In fact, a ribosome is like a mobile "workshop" which 
moves on an mRNA track and provides a platform where 
a coordinated action of many devices take place for the 
synthesis of each of the proteins. 

Only a few papers over the last few years have re- 
ported results of single-ribosome imaging and manipula- 
tion pj], El GJ, E El, [H . These experiments have es- 
tablished that in each mechano-chemical cycle, the dwell 
time of a ribosome at any codon is random. Moreover, 
this dwell time is a sum of two time intervals, namely, (i) 
the duration for which it makes a mechanical pause and 
(ii) the time it takes to translocate to the next codon. 

In this letter we report our theoretical results on 
the dwell time distribution [22| . which characterizes the 
stochastic translocation-and-pause dynamics of the ri- 
bosomes. We also introduce a few new statistical dis- 
tributions which characterize some other aspects of the 
stochastic nature of translation. We compute all these 
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statistical distributions by carrying out computer sim- 
ulations of a model of protein synthesis that captures 
both the mechano-chemistry of each individual ribo- 
some, as it moves on the mRNA template, as well as 
their in-situ steric interactions [16| . To our knowl- 
edge, none of the earlier models of "ribosome-traffic" 
d M 0, H, S OS El GJ, 0, EI have been used so far 
to investigate any of these quantitative measures of the 
stochastic kinetics of ribosomes. We treat the widths of 
these distributions as quantitative measures of "trans- 
lational noise" arising from a single mRNA template. 
Moreover, in this letter, we demonstrate the effects of 
the heterogeneity of the codon sequence of real genes on 
this "translational noise" . 

We represent the single-stranded mRNA chain by a 
one-dimensional lattice where each site corresponds a sin- 
gle codon (triplet of nucleotides). The sites i — 1 and 
i = L represent the start codon and stop codon, respec- 
tively. Each ribosome covers I sites (i.e., I codons) at a 
time; no lattice site is allowed to be covered simultane- 
ously by more than one overlapping ribosome because of 
their steric exclusion. Irrespective of the length I, each 
ribosome moves forward by only one site in each step as 
it must translate successive codons one by one. We de- 
note the position of a ribosome by the integer index of 
the leftmost lattice site it covers. 

The fig Q] captures the mechano-chemical cycle of each 
ribosome in the stage of elongation of the protein. The 
arrival of the correct amino-acid (bound to an adapter 
molecule called tRNA) and its recognition by the ribo- 
some located at the site i triggers transition from the 
chemical state 1 to 2 in the same location. The transition 
from state 2 to state 3 is driven by hydrolysis of GTP. 
Departure of the phosphate group, which is one of the 
products of GTP hydrolysis, results in the intermediate 
state 4. The peptide bond formation between the grow- 
ing protein and the incoming amino acid monomer (and 
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FIG. 1: A schematic representation of the biochemical cycle of 
a single ribosome during the elongation stage of translation 
in our model [jq . Each circle labelled by an integer index 
represents a distinct state in the mechano-chemical state of 
a ribosome. The index below the box labels the codon on 
the mRNA with which the ribosome binds. The symbols ac- 
companied by the arrows define the rate constants for the 
corresponding transitions from one state to another. 
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FIG. 2: Probability distribution of the dwell times of the ri- 
bosomes for a hypothetical homogeneous mRNA template for 
t — 12. The same set of data fits equally well with difference 
of two exponentials and with a gamma distribution (shown 
in the inset) Inset of the inset shows the coverage density 
profile of the ribosomes on the mRNA. The parameters are 
io a = uj g = 25 s" 1 , a = 0.0001 s" 1 . 



some associated biochemical processes), which leads to 
the elongation of the protein by one amino acid monomer, 
is captured by the next transition to the state 5. All 
the subsequent processes, including hydrolysis of another 
GTP molecule, the forward translocation of the ribosome 
by one codon and the departure of a naked tRNA from 
the ribosome complex are captured by a single effective 
transition from state 5 at site i to the state 1 at the site 
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FIG. 3: Typical time series of the translation events for (a) err 
gene of Escherichia coli K-12 strain MG1655 and (b) the cor- 
responding hypothetical homogeneous mRNA template, both 
corresponding to u a = 2.5s" 1 , io g — 2.5s" 1 , lo^ = 10s _1 and 
a = 0.1 s _1 . 



i + 1. More detailed explanations of the states and the 
transitions are given in ref.[16l|. 

The average number of ribosomes crossing the stop 
codon, per unit time, on the template mRNA is called 
the flux of ribosomes. The average rate of elongation of 
a protein is proportional to the average velocity of a ri- 
bosome and, therefore, the flux is a measure of the total 
rate of synthesis of the protein encoded by the mRNA 
on which the ribosomes move. The flux and the average 
density profiles of the ribosomes on the mRNA track in 
our model have been reported in ref. 16] . 

The time interval td between the arrival of a ribosome 
at a specific codon and its subsequent departure from 
there is defined as the dwell time at that codon. The run 
time T of a ribosome is the time is takes to run from the 
start codon to the stop codon on the mRNA. In other 
words, T is the time taken by a ribosome to synthesize 
a single protein. Similarly, following the terminology of 
traffic science [23| . we identify the time interval between 
the departure of the successive ribosomes from the stop 
codon as the time-headway r. Equivalently, r is the time 
interval in between the completion of the synthesis of 
successive proteins from the same mRNA template. 

In this letter we compute the distributions P(td), 
P(T), and V{t) of the probabilities of t d , T and t. We 
treat the fluctuations, i.e., root-mean-square (rms) devi- 
ations, of td, T and t as quantitative measures of noise in 
the translation of a single mRNA. Analogous measures 
of transcriptional noise have been introduced recently to 
characterize the stochasticity of polymerization of RNA 
molecules from a DNA template [26j | . 

All the calculations reported in this paper have been 
obtained by imposing open boundary conditions which 
mimics protein synthesis more realistically. The symbols 
a and f3 denote the probabilities of attachment and de- 
tachment, respectively, in time At. So, the probability 
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FIG. 4: Probability distribution of the times taken to com- 
plete the synthesis of a single polypeptide (which is identical 
to the probability distribution of the run times of ribosomes) 
for (a) err gene of Escherichia coli K-12 strain MG1655, and 
(b) the corresponding hypothetical homogeneous mRNA tem- 
plate. Both in (a) and (b), different curves correspond to 
different values of u a , all for i — 12. The discrete data 
points have been obtained from our computer simulations of 
the model whereas the lines denote the gaussian best fits to 
these data. The insets show the exponential decrease of the 
corresponding noise strengths with uj a . In both (a) and (b), 
uj g — 2.5 s _1 and a = 0.1 s -1 . 



FIG. 5: Probability distribution of the time gaps between 
the completions of the synthesis of a successive polypeptides 
(which is identical to the probability distribution of the time- 
headways in the ribosome traffic) for (a) err gene of Es- 
cherichia coli K-12 strain MG1655, and (b) the corresponding 
hypothetical homogeneous mRNA template. Both in (a) and 
(b) different curves correspond to different values of w a , all for 
I — 12. The discrete data points have been obtained from our 
computer simulations of the model whereas the lines denote 
the gamma distributions fitted to these data. The insets show 
the exponential decrease of the corresponding noise strengths 
with uj a - In both (a) and (b), ui g — 2.5 s" 1 and a = 0.1 s _1 . 



of attachment per unit time (which we call u> a ) is the 
solution of the equation a = 1 — e - w £» xAt (i n all our nu- 
merical calculations we take At = 0.001 s). Similarly, we 
define the corresponding parameter ujp for termination. 
For the same reasons as elaborated in ref . [l6| , we assume 
that ujhi — = Wh- Moreover, throughout this letter 
we use ujh = 10 s _1 , uj p = 0.0028 s" 1 , ki — 2.4 s _1 and 
0=1 s _1 which were used in ref. [la] for the bacteria 
E-coli; the values of the other parameters will be given 
in the appropriate figure captions. We incorporate the 
effects of the inhomogeneity of the sequence of codons 
in the err gene of Escherichia coli K-12 strain MG1655 
[27l | in our model exactly the same way as it was done in 



ref.[3. 

Because of the intrinsic stochasticity of the steps of the 
mechano-chemical cycle of the ribosome, the dwell time 
fluctuates even if all the codons on the mRNA track are 
identical. A typical distribution of the dwell times of 
the ribosomes during the translation of a hypothetical 
homogeneous mRNA is shown in figE] The numerical 
data obtained from computer simulations of our model 
can be fitted to the difference of two exponentials; this 
is consistent with the corresponding recent experimen- 
tal observation [22]. However, the same data fits with a 
gamma distribution almost equally well (see the inset of 
fig©- 

Typical time series of the translation events is shown 
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in fig[3] for the err gene of Escherichia coli K-12 strain 
MG1655 together with a time-series for the correspond- 
ing homogeneous mRNA template where all the rate con- 
stants other than u> a are same. The longer gaps between 
the events for the real gene arises from the fact that a 
ribosome has to wait for long periods at the "hungry 
codons" 0. 

We have plotted the distribution P(T) for the err gene 
of the Escherichia coli K-12 strain MG1655, for differ- 
ent values of the model parameters cu a in figHl the data 
for the corresponding hypothetical homogeneous mRNA 
template are plotted in figQJb). In figJ5] we have plot- 
ted the corresponding data for V T . The variation of 
the strength of the noise with the model parameters are 
shown in the insets of the respective figures. Both the 
measures of translational noise fall exponentially with the 
increase of uj a . In other words, increase in the availabil- 
ity of the monomeric subunits (which is indicated by to a ) 
reduce the noise level. Similar trend of variation of noise 
with u>h (i.e., the rate of "fuel" consumption) has been 
observed (but not shown graphically). 

Comparing the data in figllja) and fig[4jb) we conclude 
that the sequence inhomogeneity of real genes not only 
slows down the polymerization of the proteins, but also 
makes the process more noisy as compared to the trans- 
lation of the hypothetical homogeneous gene. Similarly, 
comparing the data in figEJa) with those in figlDJb) we 
establish that sequence inhomogeneity of real genes leads 
to longer mean, as well as stronger fluctuations, in r than 



for the hypothetical homogeneous template. 

The data for P(T), obtained from computer simula- 
tions, fit well with a gaussian distribution In contrast, 
the best fit to those for V T is a gamma disribution Such 
long-tail distributions are quite common in gene expres- 
sion and describe the characteristic features of various 
statistical properties of gene expression [H, [29|, [3(| ■ 

In this letter we have developed a new conceptual 
framework for analyzing the intrinsic stochasticity in the 
process of polymerization of proteins by ribosome ma- 
chines from a single mRNA template. The widths of 
the statistical distributions, which characterize different 
aspects of this stochasticity, serve as quantitative mea- 
sures of noise in the translation of a single mRNA. By 
comparing our results for a specific gene of the bacte- 
ria Escherichia coli with those for the corresponding ar- 
tificial homogeneous mRNA template, we have demon- 
strated the effects of the sequence inhomogeneities of real 
genes on the translational noise. The nature of the dwell 
time distributions predicted by our theory is consistent 
with the corresponding observations [22j in recent single- 
ribosome experiments. We hope our other predictions 
will stimulate new experiments on translational noise. 
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